Quantum Critical Transport Near the Mott Transition 
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We perform a systematic study of incoherent transport in the high temperature crossover region 
of the half-filled one-band Hubbard model. We demonstrate that the family of resistivity curves 
displays characteristic quantum critical scaling of the form p(T,6U) — p c (T)f(T/T (SU)), with 
T (5U) \SU\ ZU , and p c (T) ~ T. The corresponding /3-function displays a "strong coupling" form 
P ~ ln(p c /p), reflecting the peculiar mirror symmetry of the scaling curves. This behavior, which is 
surprisingly similar to some experimental findings, indicates that Mott quantum criticality may be 
acting as the fundamental mechanism behind the unusual transport phenomena in many systems 
near the metal-insulator transition. 
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Many systems close to the metal-insulator transition 
(MIT) often display surprisingly similar transport fea- 
tures in the high temperature regime [THS]. Here, the 
family of resistivity curves typically assumes a char- 
acteristic "fan-shaped" form (see Fig. 1(a)), reflecting 
a gradual crossover from metallic to insulating trans- 
port. At the highest temperatures the resistivity de- 
pends only weakly on the control parameter (concentra- 
tion of charge carriers pQ, or pressure [2j |3]), while as 
T is lowered, the system seems to "make up its mind" 
and rapidly converges towards either a metallic or an in- 
sulating state. Since temperature acts as a natural cut- 
off scale for the metal-insulator transition, such behav- 
ior is precisely what one expects for quantum criticality 
[4]. In some cases [I], the entire family of curves dis- 
plays beautiful scaling behavior, with a remarkable "mir- 
ror symmetry" of the relevant scaling functions [4]. But 
under which microscopic conditions should one expect 
such scaling phenomenology? What is the corresponding 
driving force for the transitions? Despite recent progress, 
such basic physics questions remain the subject of much 
ongoing controversy and debate. 

The phenomenon of disordered-driven Anderson local- 
ization of noninteracting electrons is at present rather 
well understood based on the scaling formulation [5 J , and 
is generally viewed as an example of a T = quantum 
phase transition (QPT). On the other hand, a consid- 
erable number of recent experiments [lj provide com- 
pelling evidence that strong correlation effects - some 
form of Mott localization - may be the dominant mech- 
anism |6J. Should one expect similar or very different 
transport phenomenology in the Mott picture? Is the 
paradigm of quantum criticality even a useful language 
to describe high temperature transport around the Mott 
point? These issues are notoriously difficult to address, 
because conventional Fermi liquid concepts simply can- 
not be utilized in the relevant high temperature incoher- 
ent regime. In this Letter we answer this question in the 




Figure 1: (Color online) (a) DMFT resistivity curves as func- 
tion of temperature along different trajectories —0.2 < SU < 
+0.2 with respect to the instability line SU = (black dashed 
line, see the text). Data are obtained using IPT impurity 
solver, (b) Resistivity scaling; essentially identical scaling 
functions are found from CTQMC (open symbols) and from 
IPT (closed symbols) 

framework of dynamical mean- field theory (DMFT) [1], 
the only theoretical method that is most reliable precisely 
at high temperatures. 

Model and DMFT solution. We consider a single-band 
Hubbard model at half-filling 
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where c\ a and Ci a are the electron creation and annihi- 
lation operators, = c\ a Ci a , tij is the hopping ampli- 
tude, and U is the repulsion between two electrons on 
the same site. We use a semicircular density of states, 
and the corresponding half-bandwidth D is set to be our 
energy unit. We focus on the paramagnetic DMFT so- 
lution, which is formally exact in the limit of large co- 
ordination. Here the Hubbard model maps onto an ef- 



2 



i i i | i i i i | i i i i | y i i | i 1 i i | i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i 




Figure 2: (Color online) DMFT phase diagram of the fully 
frustrated half-filled Hubbard model, with a shaded region 
showing where quantum critical- like scaling is found. Metal- 
lic U C 2(T) and insulating U c i(T) spinodals (dotted lines) are 
found at T < T c ; the corresponding first order phase tran- 
sition is shown by a thick solid line. The thick dashed line, 
which extends at T > T c shows the "instability trajectory" 
U* (T) , and the crossover temperature T Q delimits the QC re- 
gion (dash-dotted lines). The inset shows examples of eigen- 
value curves at three different temperatures, with pronounced 
minima at U*(T) determining the instability trajectory. 

fective Anderson impurity model supplemented by a self- 
consistency condition [lj. To solve the DMFT equations 
we use the iterated perturbation theory (IPT) [lj and 
cross-check our results with numerically exact continuous 
time quantum Monte Carlo (CTQMC) 019]. We find, in 
agreement with previous work [3j, that after appropriate 
energy rescaling (see below), the two methods produce 
qualitatively and even quantitatively identical results in 
the incoherent crossover region that we examine. 

It is well known that at very low temperatures T < 
T c ~ 0.03, this model features a first order metal- 
insulator transition terminating at the critical end-point 
T c (Fig. 2), very similar to the familiar liquid-gas transi- 
tion [3]. For T > T c , however, different crossover regimes 
have been tentatively identified [TJ [11] but they have not 
been studied in any appreciable detail. The fact that the 
first order coexistence region is restricted to such very 
low temperatures provides strong motivation to examine 
the high temperature crossover region from the perspec- 
tive of "hidden quantum criticality". In other words, the 
presence of a coexistence dome at T < T c <C 1, an effect 
with very small energy scale, is not likely to influence 
the behavior at much higher temperatures T ^> T c . In 
this high temperature regime smooth crossover is found, 
which may display behavior consistent with the presence 
of a "hidden" quantum critical (QC) point at T = 0. To 
test this idea, we utilize standard scaling methods appro- 
priate for quantum criticality, and compute the resistiv- 
ity curves along judiciously chosen trajectories respecting 



the symmetries of the problem. 

Instability trajectory formalism. Previous work has al- 
ready recognized [3] that, in order to reveal the proper 
scaling behavior close to the critical end-point, one has to 
follow a set of trajectories parallel to "zero field" trajec- 
tory U*(T). We thus expect SU = U-U*(T) to play the 
role of the scaling variable corresponding to a symmetry- 
breaking field favoring one of the two competing (metal 
vs. insulator) phases. By analogy [3j H2] to the familiar 
liquid-gas transition, we determine the precise location 
of such an "instability trajectory" by examining the cur- 
vature of the corresponding free energy functional [T3] . 
This curvature vanishes at T c and is finite and minimal 
at T > T c , along this instability line. Consequently, as 
in Refs. [3j [131 EL our problem is recast as an eigen- 
value analysis of the corresponding free energy functional 
F[G(iw n )] for which the DMFT Green's function solu- 
tion G jj m ft{^ w n) represents a local extremum, and can 
be regarded as a vector in an appropriate Hilbert space. 

The free energy near such an extremum can be written 
as F[G(iw n )\ = To + Tt 2 £ m , n 5G(iu m )M mn 5G(iu> n ) + 
.... where 
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and 5G(iuj n ) = G(iuj n ) — GoMFT^n)- The curvature 
of the free energy functional is determined by the lowest 
eigenvalue A of the fluctuation matrix M. As explained 
in Supplementary Notes [15J, A can be obtained from the 
iterative solution of DMFT equations. The difference of 
the Green's functions in iterations n and n + 1 of the 
DMFT self-consistency loop is given by 

£G (n+1) (iuj n ) - SG M (iuj n ) = e- nX SG^ (iuj n ) (3) 

and therefore A determines the rate of convergence of the 
Green function to its solution. 

An example of our calculations is shown in the inset of 
Fig. 2, where the eigenvalues at several temperatures are 
plotted as a function of interaction U/U c . The minima 
of these curves define the locus of the " instability trajec- 
tory" U*(T) , which terminates at the critical end-point 
(J7 C ,T C ), as shown on Fig. 2. Note, that the immediate 
vicinity T w T c of the critical end-point has been care- 
fully studied theoretically [3j, and even observed in ex- 
periments [2j, revealing classical Ising scaling (since one 
has a finite temperature critical point) of transport in this 
regime. In our study, we examine the crossover behav- 
ior at much higher temperatures T ^> T c , displaying very 
different behavior: precisely what is expected in presence 
of quantum criticality. 

Resistivity calculation. To reveal quantum critical scal- 
ing, we calculate the temperature dependence of the resis- 
tivity along a set of trajectories parallel to our instability 
trajectory (fixed 5U = U — U*(T)). Resistivity was cal- 
culated using standard DMFT procedures [l], with the 
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Figure 3: (Color online) (a) Scaling parameter To(5U) as func- 
tion of control parameter SU = U — U*\ the inset illustrates 
power-law dependence of scaling parameter T = c\SU\ zu ] (b) 
Resistivity p c (T) of the "separatrix". Excellent agreement is 
found between IPT (closed symbols) and CTQMC (open sym- 
bols) results. 

maximum entropy method [16J utilized to analytically 
continue CTQMC data to the real axis. The resistivity 
results are shown in Fig. 1, where on panel (a) IPT resis- 
tivity data for SU = 0, ±0.025, ±0.05, ±0.1, ±0.15, ±0.2 
in the temperature range of T « 0.07 — 0.2 are presented 
(CTQMC data are not shown for the sake of clarity of the 
figure). The resistivity is given in units of pMotu maxi- 
mal resistivity according to the Boltzmann quasi-classical 
theory of transport [IT]- The family of resistivity curves 
above (SU > 0) the "separatrix" p c (T) (dashed line, cor- 
responding to SU = 0) has an insulating-like behavior, 
while metallic dependence is obtained for SU < 0. 

Scaling analysis. According to what is generally ex- 
pected for quantum criticality, our family of curves 
should satisfy the following scaling relation: 



p(T,SU) = p c (T)f(T/T (SU)). 



(4) 



We thus first divide each resistivity curve by the "separa- 
trix" p c (T) = p(T, SU = 0) and then rescale the temper- 
ature, for each curve, with an appropriately chosen pa- 
rameter To (SU) to collapse our data onto two branches, 
Fig. 1(b). Note that this unbiased analysis does not 
assume any specific form of T Q (SU): it is determined for 
each curve simply to obtain optimum collapse of the data 
[18J. This puts us in a position to perform a stringent test 
of our scaling hypothesis: true quantum criticality corre- 
sponds to T Q (SU) which vanishes at SU = and displays 
power law scaling with the same exponents for both scal- 
ing branches. As seen in Fig. 3(a) To falls sharply as 
U = U* is approached, consistent with the QC scenario, 
but opposite to what is expected in a "classical" phase 
transition. The inset of Fig. 3(a) with log- log scale shows 
clearly a power-law behavior of T = c\SU\ zu , with the es- 

IPT 

timated power (zv) = 0.56 ±0.01 for "metallic " side, 

IPT 

and (zv) = 0.57 ± 0.01 for an insulating branch. 
We also find (Fig. 3(b)) a very unusual form of our 



critical resistivity p c (T), corresponding to the instability 
trajectory. Its values largely exceeds the Mott limit, yet 
it displays metallic-like but non-Fermi liquid-like tem- 
perature dependence p c (T) ~ T. Such puzzling behavior, 
while inconsistent with any conventional transport mech- 
anism, has been observed in several strongly correlated 
materials close to the Mott transition [T71 [20] . Our re- 
sults thus suggest that it represents a generic feature of 
Mott quantum criticality. 

f3— function and mirror symmetry of scaled curves. To 
specify the scaling behavior even more precisely, we com- 
pute the corresponding /3-function [4] /3(g) = with 
g = pcj p being the inverse resistivity scaling function. 
Remarkably (Fig. 4(a)), it displays a nearly linear de- 
pendence on Ing, and is continuous through SU = in- 
dicating precisely the same form of the scaling function 
on both sides of the transition - another feature exactly of 
the form expected for genuine quantum criticality. This 
functional form is very natural for the insulating trans- 
port, as it is obtained even for simple activated behavior 
p(T) ~ e~ E ^ T . The fact that the same functional form 
persists well into the metallic side is a surprise, especially 
since it covers almost an order of magnitude for the resis- 
tivity ratio. Such a behavior has been interpreted |4J to 
reflect the "strong coupling" nature of the critical point, 
which presumably is governed by the same physical pro- 
cesses that dominate the insulator. This points to the 
fact that our QC behavior has a strong coupling, i.e. 
non- perturb at ive character. 

The fact that the j3— function assumes this logarithmic 
form on both sides of the transition is mathematically 
equivalent |4J to stating that the two branches of the 
corresponding scaling functions display "mirror symme- 
try" over the same resistivity range. Indeed, we find that 
transport in this QC region exhibits a surprisingly de- 
veloped reflection symmetry (dash vertical lines of Fig. 
4(a) mark its boundaries). Such a symmetry is clearly 
seen in Fig. 4(b), where the resistivity p/p c (for SU > 0) 
and conductivity a/a c = p c /p (SU < 0) can be mapped 



onto each other by reflection with p ^ u ^ 



aj-SU) 
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Note that T/T Q — 1 sets the boundary of the quan- 
tum critical region, over which the reflection symmetry of 
scaled curves is observed. It is depicted by dash-dotted 
crossover lines To on phase diagram of Fig. 2 [15] . 

These remarkable features of the f3— function, and asso- 
ciated reflection-symmetry, have been observed earlier in 
experimental [TJ|2T] and theoretical [4] studies, which ten- 
tatively associated this with disorder-dominated MITs. 
Speculation that f3 ~ In g reveals disorder as the funda- 
mental driving force for MIT, presumably reflects the fact 
that, historically, it has first been recognized for Ander- 
son transitions [5j. Our work, however, shows that such 
behavior can be found even in absence of disorder - in in- 
teraction driven MITs. This finding calls for re-thinking 
of basic physical processes that can drive the MIT. 



4 




ln(p c /p) T/T Q 

Figure 4: (Color online) (a) f3— function shows linear in 
ln(p c /p) behavior close to the transition. Open symbols are 
for metallic branch (5U < 0) and closed ones are for the insu- 
lating side (SU > 0) ; vertical dashed lines indicate the region 
where mirror symmetry of curves is found, (b) Reflection 
symmetry of scaled curves close to the transition. 

Conclusions: We have presented a careful and de- 
tailed study of incoherent transport in the high tempera- 
ture crossover regime above the critical endpoint T c of a 
single-band Hubbard model. Our analysis revealed a so- 
far overlooked scaling behavior of the resistivity curves, 
which we interpreted as evidence of "hidden" Mott quan- 
tum criticality. Precisely locating the proposed QCP in 
our model is hindered by presence of the low tempera- 
ture coexistence dome, which limits our quantum critical 
scaling to the region well above T c . Regarding the nature 
of transport in the QC regime, we found that the critical 
resistivity well exceeds the Mott limit, and yet it - sur- 
prisingly - assumes a metallic form, in dramatic contrast 
to conventional MIT scenarios. These features, together 
with large amounts of entropy characterizing this entire 
regime [22J, prove surprisingly reminiscent of the "holo- 
graphic duality" scenario [23| [24] for a yet-unspecified 
QC point. Interestingly, the holographic duality picture 
has - so far - been discussed mostly in the context of 
quantum criticality in correlated metals (e.g. T = mag- 
netic transitions in heavy fermion compounds). Ours is 
the first work proposing that the same physical picture 
could apply to quantum criticality found at the MIT. 

We believe that our results provide a significant new 
perspective on QC around the Mott transition, and a 
deeper understanding of an apparent universality in the 
high temperature crossover regime. Our method traces 
a clear avenue for further searches for QC scaling, which 
are likely to be found in many other regimes and models. 

In particular, it would be interesting to study a cor- 
responding critical regime by going beyond the single- 
site DMFT analysis. It was shown in [19J that inclusion 
of spatial fluctuations does not significantly modify the 
high temperature crossover region in half-filled Hubbard 
model. Consequently, we expect our main findings to 
persist. An even more stringent test of our ideas should 
be provided in models where the critical end-point T c 



can be significantly reduced. This may include studies of 
the Mott transition away from half-fillling [25J, or in sys- 
tems with frustrations [6j [26] . In such situations the pro- 
posed scaling regime should extend to much lower tem- 
peratures, perhaps revealing more direct evidence of the 
- so far - hidden Mott QC point. Our ideas should also 
be tested by performing more detailed transport experi- 
ments in the relevant incoherent regime, a task that may 
be easily accessible in various organic Mott systems [3], 
where T c is sufficiently below room temperature. 
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SUPPLEMENTARY NOTES: QUANTUM 
CRITICAL TRANSPORT NEAR THE MOTT 
TRANSITION 

Eigenvalue analysis of the free energy curvature 

Here we present in detail the procedure that we use 
to determine the minimum curvature of the free energy 
functional for a given temperature. For simplicity we 
concentrate on the Bethe lattice. The Ginzburg-Landau 
free energy functional F{G) in the Hilbert space of the 
Matsubara Green's functions G = G(iuo n ) is given by [1- 

m 



F{G) = —Tt 2 G 2 + F imp (G) 

= -Tt 2 Y,G 2 (iuj n ) + F imp (G), 



(5) 



where the first term is the energy cost of forming the 
Weiss field A = t 2 G around a given site, while the second 
term is the free energy of an electron at this site in the 
presence of the Weiss field. 

Close to a local extremum Go, we can Taylor expand 
F(G) in terms of deviation from this point 5G = G — Go: 



F(G) = F(G ) + Tt 2 Y,SG(ioj rn )M mn 5G(iuj n ) ■ 



F(G )+T£ 2 5GM5G + ... 



where 



1 



d 2 F(G) 



2Tt 2 dG(iuj n )dG(iuj m ) 
We introduce a gradient function 



G=G 



9(G) 



1 dF(G) 



MSG 



2Tt 2 dG 

and define an iteration-substitution procedure by 
g&n+i) =SG M -g(G M ), 



(6) 



(7) 



(8) 



(9) 



which gives the minimum of the free energy as the itera- 
tion step n — > oo. In the case of the free energy functional 
given by Eq. (J5|, we have 



g(G)=G imp (G)-G 



(10) 



and in the iterative procedure the Green's function con- 
verges to the minimum of the free energy which is at the 
same time also a self-consistent solution of the DMFT 
equations, given by the relation Gi mp (iuj n ) = G(iuj n ). 

The curvature of the free energy for interaction U and 
temperature T can be obtained as follows. The eigenbasis 
G a and eigenvalues X a of matrix M are defined by 
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Figure 5: Convergence rate in the iterative solution of DMFT 
equations at T/T c = 2.33 using IPT impurity solver, panel 
(a), and CTQMC impurity solver, panel (b). The dashed 
lines in panel (b) are linear fits to the data. The insets are 
the corresponding eigenvalues determined by the slopes from 
the main panels. 



We can expand SG^ as 



(12) 



where are the coefficients of SG^ in the eigenvalue 
basis. Substituting into Eq. (l9|), one obtains 



where 



ln(l-A Q ). 



(13) 



(14) 



For large n the term with lowest B a = B ao , which cor- 
responds to the lowest eigenvalue X ao = A, is dominant 



SG M 



)aL 0) G, 



a ^ao ' 



n > 1. 



(15) 



MG a — X a G a . 



(ii) 



Here is the coefficient corresponding to the Green 
function with the lowest eigenvalue A. Now it is obvious 
that through iterations, the solution G approaches to Go 
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Figure 6: Phase diagram obtained with IPT, panel (a), and 
CTQMC, panel (b). Temperature and interaction are scaled 



by their values at the critical endpoint (T c 



0.046, 



rpQMC = 03 and jjIpt = 2 4725 jjQmc = 2 Red dashed 

line is the instability line U* (T) , full red line is the line of the 
first order MIT, and green and blue dotted lines are left and 
right spinodals. 



exponentially along a direction defined by the eigenvec- 
tor of M corresponding to its minimal eigenvalue A. The 
coefficient B ao and the corresponding eigenvalue A are 
then obtained from the slope in the iterative relation 



In 



G(iuj n )^ - G(iu n ) 



(n) _ 



const — nB a 



(16) 



which follows from Eq. ([13]). 

In practice, to obtain A (and thus the curvature of 
free energy), we monitor DMFT loop convergence rate, 
G(icj )( n+1 ) — G(iuj )^ n \ in as many iterations as possible 
and then linearly fit In {G(iuj )( n+1 ^ - G(iu )^) versus 
iteration index n. Here uo Q — ttT is the lowest Matsubara 
frequency. For small A, B ao « A. We repeat this proce- 
dure for different values of U at the same temperature T 
to determine U*(T) in which X(U)\t is minimal. It takes 
few iterations of the DMFT loop to enter into the linear 



regime given by Eq. (16). 

With IPT impurity solver, we can use data from several 
tens of iterations to determine the slope B ao , Fig. 1(a). 
The solution with CTQMC impurity solver has a statisti- 
cal error and the number of iterations is limited before the 
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Figure 7: (a) Resistivity as function of temperature for 
SU = ±0.025,0.05,0.1,0.15,0.2. (b) By scaling the data along 
T— axis by To data are collapsed onto two branches. Data are 
collapsed on the lowest SU = ±0.025 curves. 



difference |G(zco> )( n+1 ) — G(iuj )^\ becomes too small 
and acquires a large relative statistical error. Neverthe- 
less, we were able to determine rather precisely the eigen- 
value A and the interaction U*(T) for which it becomes 
minimal, Fig. 1(b). The "instability line" corresponding 
to the minimum curvature of the free energy is shown in 
Fig. 2(a) (IPT phase diagram), and Fig. 2(b) (CTQMC 
phase diagram). Error bars in Fig. 2(b) are estimates of 
the uncertainty in the position of the instability line. 



Details of the scaling procedure 

The resistivity p(T, SU) is calculated along the lines 
parallel to the instability line U*(T). Here SU = U - 
U*(T). The resistivity is fist divided by its value p c (T) at 
SU = 0, Fig. 3(a). Then for each SU the temperature axis 
is scaled by To where the scaling parameter To is chosen 
to collapse the data onto two branches: insulating-like 
for SU > and metallic-like for SU < 0. The scaling was 
done in such a way that data were collapsed on the lowest 
curves with SU = ±0.025 as shown in Fig. 3(b). The 
scaling parameter To has a power law form To = c\SU\ zv , 
where the prefactor c depends on this referent value of 
U. 

Our data exhibit a reflection symmetry which is seen 
in Fig. 4(a), where the resistivity p/p c (for SU > 0) and 
conductivity cr/a c = p c / p (SU < 0) can be mapped onto 
each other by reflection with p ( su ^ = Pc(-su) ^ 

The standard estimate for the scale (prefactor) of the 
crossover temperature is obtained by requiring that the 
scaling variable x = T/To = 1 at the point where 
the scaling function changes its functional form; in our 
case this corresponds to the temperature below which 
the mirror symmetry of the scaling curves no longer 
holds. Before rescaling the prefactor, this is found at 
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Figure 8: (a) Total scaled data (for clarity shown by single 
branches) exhibit reflection symmetry for T > T sc . (b) T/Tq- 
axis is rescaled by T sc in such a way that T/T Q = 1 set the 
boundary of quantum critical region over which the reflection 
symmetry of scaled curves is observed. 
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Figure 9: (a) Scaling temperature To vs. SU obtained from 
scaling procedure shown on Fig. 4(a). (b) To vs. SU after 
rescaling by T sc . The boundary of the quantum critical region 
is now given by condition T/T G = 1. 



x* = T SC /T Q = 0.052. Our final form of scaled data 
is shown in Fig. 4(b), where T/Tq - axis is rescaled by 
x* = T SC /T Q so that T/T Q = 1 sets the boundary of the 
quantum critical region. The scaling parameter To as 
function of SU is shown in Fig. 5(a) (before rescaling 
with x*) and in Fig. 5(b) (after rescaling with x*). The 
corresponding values for c and zv from the power law fit 
are also given. 

We can now plot the crossover temperature To setting 
the boundary of QC region on our phase diagram. T = Tq 



condition is equivalent to 

T = c\8U(T o )r 



c5\U-U*(T D )r, 



(17) 



where U*(T Q ) is value of U at temperature T = T , along 
the instability line. This equation implicitly defines the 
crossover line T (U). Alternatively, we can invert this 
dependence to describe the same crossover line as Uq(T) 
that takes the form 



U±(T) = U*(T)±(T/c) 



l/zv 



(18) 



As we can see from this expression, the crossover line ap- 
proaches the instability line at low T and diverges away 
from it at high T. The phase diagram including the in- 
stability line and the crossover temperature To is shown 
on Fig. 6. 
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Figure 10: (Color online) DMFT phase diagram of the fully 
frustrated half- filled Hubbard model. The thick dashed line, 
which extends at T > T c shows the "instability trajectory" 
U*(T), and the crossover temperature T G delimits the QC 
region (dash-dotted lines). 
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